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Abstract 

Halo nuclei are a promising new arena for studies based on effective field theory (EFT). We 
develop an EFT for shallow p-w&ve states and discuss the application to elastic na scattering. In 
contrast to the s-wave case, both the scattering length and effective range enter at leading order. 
We also discuss the prospects of using EFT in the description of other halos, such as the three-body 
halo nucleus 6 He. 
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I. INTRODUCTION 



Nuclear halo states have been found in a number of light nuclei close to the nucleon drip 
lines. They are characterized by a very low separation energy of the valence nucleon (or 
cluster of nucleons). As a consequence, the nuclear radius is very large compared to the size 
of the tightly bound core. The large size of halo nuclei leads to threshold phenomena with 
important general consequences for low-energy reaction rates in nuclear astrophysics. One 
example is the reaction p + 7 Be — > 8 B + 7 which is important for solar neutrino production. 
The nucleus 8 B is believed to be a two-body proton halo, consisting of a 7 Be core and a 
proton PJ. Somewhat more complicated are three-body halos consisting of a core and two 
slightly bound nucleons. Particularly interesting are Borromean three-body halos, where no 
two-body subsystem is bound. Typical examples are 6 He and n Li, which consist of a 4 He 
and 9 Li core, respectively, and two neutrons. For reviews of halo nuclei, see Ref. 0. The 
physics of halo nuclei is an important part of the physics program at RIA Jl . A thorough 
discussion of reactions with rare isotopes can be found in Ref. [Q . 

The physics of halo nuclei is a promising new arena for Effective Field Theory (EFT). 
EFTs provide a powerful framework to explore separation of scales in physical systems in 
order to perform systematic, model- independent calculations ||. If, for example, the relative 
momentum k of two particles is much smaller than the inverse range of their interaction 1/-R, 
observables can be expanded in powers of kR. All short-distance effects are systematically 
absorbed into a few low-energy constants using renormalization. The EFT approach allows 
for systematically improvable calculations of low-energy processes with well-defined error 
estimates. The long-distance physics is included explicitly, while the corrections from short- 
distance physics are calculated in an expansion in the ratio of these two scales 1 . The inherent 
separation of length scales in halo nuclei makes them an ideal playing ground for EFT. 

In recent years, there has been much interest in applying EFT methods to nuclear systems 
M,M- Up to now, nuclear EFT has mainly been applied to two-, three-, and four-nucleon sys- 
tems starting from a fundamental nucleon-nucleon interaction. The original motivation was 
to understand the gross features of nuclear systems from a QCD perspective by deriving the 
nuclear potential and currents relevant for momenta comparable to the pion mass (p ~ m n ) 
||. More recently, it has been realized that it is possible to carry out very precise calcula- 
tions for fundamental physics processes at lower energies. For very low momenta (p^m^), 
even pion exchange can be considered "short-distance" physics. In this case, one can use an 
effective Lagrangian including only contact interactions. The large s-wave scattering lengths 
require that the leading two-body contact interaction be treated nonperturbatively || |10 |. 



In the two-nucleon system, this program has been very successful (see, e.g., Refs. |6|, |TT| and 
references therein). 

Using EFT, one can relate low-energy measurements in one reaction to observables in a 
similar (but unmeasured) reaction in a controlled expansion with reliable error estimates. 
This is in contrast to standard potential model calculations where errors can only be esti- 
mated by comparing different potentials. An example of a precise calculation in the "pion- 
less" EFT is the reaction n + p — + d + 7 [Q, which is relevant to big-bang nucleosynthesis 



Note that "effective theory" is sometimes used in reference to a model that captures the essence of the 
relevant long-distance physics without necessarily accounting for the short-distance physics in a systematic 
way. Here we use "EFT" in the model-independent sense described above, in which "power counting" of 
the different orders in the expansion is a crucial ingredient. 
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(BBN). As for many other reactions of astrophysical interest, the uncertainty in the cross 
section is difficult to determine due to the lack of data at low energies and the lack of in- 
formation about theoretical estimates. In the energy of relevance to BBN, both El and 
Ml capture are important. They have been calculated to fifth and third order, respectively, 
where two new counterterms appear. Using the measured cold-capture rate and data for 
the deuteron photodisintegration reaction to fix the counterterms, the n + p — > d + 7 cross 
section was computed to 1% for center-of-mass energies E < 1 MeV. 

Much of the strength of EFT lies in the fact that it can be applied without off-shell 
ambiguities to systems with more nucleons. The crucial issue of the relative size of three- 



body forces has been investigated in the three-body system fl3fl . Nucleon-deuteron scattering 



in all channels except the S1/2 wave can be calculated to high orders using two-nucleon input 



only, with results in striking agreement with data and potential-model calculations H14] . For 
example, to third order, the s 3 / 2 scattering length is found to be a^f 7 ^ = 6.33 ± 0.10 fm, 

to be compared to the measured afe = ^-35 ± 0.02 fm. In contrast, in the Si/ 2 -wave 
channel, the non-perturbative running of the renormalization group requires a momentum- 
independent three-body force in leading order ||15|| . Once the new parameter is fitted to 
(say) the scattering length, the energy dependence is predicted. The triton binding energy, 
for example, is found to be = 8.0 MeV in leading order, already pretty close to the 

experimental B^ xp€) = 8.5 MeV. Recently, this approach was also applied to Ad scattering 



and the hypertriton [|TjJ. Using the hypertriton binding energy to fix the three-body force, 
the low-energy Ad scattering observables can be predicted. The results are very insensitive 
to the poorly known AN low-energy parameters. In a related study, the Ad Phillips line was 
established |~7| . 

However, in an EFT it is by no means necessary to start from a fundamental nucleon- 
nucleon interaction. If, as in halo nuclei, the core is much more tightly bound than the 
remaining nucleons, it can be treated as an explicit degree of freedom. One can write an 
EFT for the contact interactions of the nucleons with the core and include the substructure of 
the core perturbatively in a controlled expansion. This approach is appropriate for energies 
smaller than the excitation energy of the core. In other words, one can account for the 
spatial extension of the core by treating it as a point particle with corrections from its 
finite size entering in a derivative expansion of the interaction. This is a consequence of the 
limited resolution of a long wavelength probe which cannot distinguish between a point and 
an extended particle of size R if the wavelength A ^> R. 

In this paper, we consider the virtual p-wave state in na scattering as a test case. Even 
though there is no bound state in this channel, it has all the characteristics of a two-body 
halo nucleus. Furthermore it is relevant for the study of the Borromean three-body halo 
6 He, which will be addressed in a forthcoming publication [JTSJ . Elastic na scattering is 
relatively well known experimentally. Since the nucleon has j = 1/2 and the a particle 
has j = 0, there are contributions from an s wave (si/2), two p waves (pi/2 and P3/2), etc.. 
Arndt, Long, and Roper performed a phase-shift analysis of low-energy data and extracted 
the effective range parameters in the s and p waves ||19|| . The p ?J / 2 partial wave displays a 
resonance at E ~ 1 MeV corresponding to a shallow virtual bound state, while the S1/2 and 
P1/2 partial waves are nonresonant at low energies. We will show that this p-wave resonance 
leads to a power counting different from the one for s-wave bound states [§, [T(J that has 



3 



been discussed extensively in the literature because of its relevance for the no-core EFT. 2 
In particular, proper renormalization requires two low-energy parameters at leading order, 
namely the scattering length and the effective range. The extension to higher-orders is 
straightforward. As we will see, the EFT describes the low-energy data very well. 

The organization of this paper is as follows: In the next Section, we work out renormaliza- 
tion and power counting for a p-wave resonance in the simpler context of spinless fermions. 
In Section |TJ, we include the spin and isospin of the nucleon and apply our formalism to 
elastic na scattering. In Section [TV], we summarize our results and present an outlook. In 
particular, we discuss the extension to the Borromean three-body halo 6 He and the reaction 
p + 7 Be -> 8 B + 7. 

II. EFT FOR SHALLOW P-WAVE STATES 

In this section, we develop the power counting for shallow p-wave states (bound states 
or virtual states) in the particularly simple case of a hypothetical system of two spinless 
fermions of common mass m. Our arguments are a generalization of those in Ref. B. 

In order to have a shallow bound state, we need at least two momentum scales: the 
breakdown scale of the EFT, M^, and a second scale, M\ <C M^, that characterizes the 
shallow bound state. The scale M hi is set by the degrees of freedom that have been integrated 
out. In the case of an EFT without explicit pions and core excitations, is the smallest 
between the pion mass m n and the momentum corresponding to the energy of the first 
excited state. The scale M\ is not a fundamental scale of the underlying theory. It can be 
understood as arising from a fine tuning of the parameters in the underlying theory. If the 
values of these parameters were changed slightly, the scale M\ would disappear. We seek 
an ordering of contributions at the scale M\ Q in powers of Mi a /Mf^. Due to the presence of 
fine-tuning, naive dimensional analysis cannot be applied. 

For simplicity we neglect relativistic corrections. They are generically small because they 
are suppressed by powers of the particle mass m, and in the cases of interest here m 3> M hi . 
They can be included along the lines detailed in Ref. ||. 

A. Natural Case 

First, we will consider the natural case without any fine-tuning. The scale of all low- 
energy parameters is then set by and naive dimensional analysis can be applied. 

The T-matrix for the non-relativistic scattering of two spinless fermions with mass m in 
the center-of-mass frame can be expanded in partial waves as 

r M = gr, ft co s9 ) ^g.^±I_P,(cose), (i) 

where k is the center-of-mass momentum, 9 the scattering angle, and Pi(cos8) is a Legendre 
polynomial. 3 The generalized effective range expansion for arbitrary angular momentum I 



2 By no-core EFT we mean an EFT where all nuclei are dynamically generated from nucleon (and possibly 
pion and delta isobar) degrees of freedom. 

3 Note that we assume the two fermions are distinguishable. 
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FIG. 1: Lowest-order diagrams for the perturbative expansion of T\. Particle propagators are 
represented by solid lines. The C\ (Cf) vertex is indicated by the open circle (crossed circle). 



reads: 



fc^cottf, 



(2) 



ai 2 4 

ri, and Vi are the scattering length, effective range, and shape parameter in 



where ai 

the l-th partial wave, respectively. For I = 0, Eq. (§) reproduces the familiar effective range 
expansion for s waves. Note that the dimension of the effective-range parameters depends on 
the partial wave. In the s wave, ao and ro have dimensions of length, while Vo has dimensions 
of (length) 3 . For p waves, a\ has dimension of (length) 3 (it is a scattering "volume"), r% has 
dimension l/(length) (it is an "effective momentum"), and V\ has dimension of length. 

The s-wave contribution T has been discussed in detail in the literature || [10). Our goal 
here is to set up an EFT that reproduces the p-wave contribution T\ in a low-momentum 
expansion, 



Ti(k, cost 



12na 



-k 2 cos6 ( 1 



rn 



^2_ miA; 3 + ^ (air 2_p i)A ,4 + 



(3) 



We start with the most general Lagrangian for spinless fermions with p-wave interactions: 
^ 2 



id + 



2m 



C p <-» <-> C p 



(^V'V^t^V^+H.c. 



+ . . . 



(4) 



where V is the Galilean invariant derivative, H.c. denotes the Hermitian conjugate, 

and the dots denote higher- derivative interactions that are suppressed at low energies. The 
fermion propagator is simply 



Po — p 2 /2m + ie 



(5) 



and the Feynman rules for the vertices can be read off Eq. (^). Around the non-relativistic 
limit, all interaction coefficients contain a common factor of 1/m that follows from Galilean 
invariance. ^From dimensional analysis, we have Cf ~ \2 / njmMf li and Cf ~ 127T y 'mM^. 
The exact relation of Cf and C\ to the scattering length and effective range will be obtained 
in the end from matching to Eq. @. 

We work in the center-of-mass frame and assign the momenta ±k and ±k' to the incoming 
and outgoing particles, respectively. The total energy is E = k 2 /m = k' 2 /m. The EFT 
expansion is in powers of k/M^. The leading contribution to T\ is of order 12irk 2 /mM^. It 
is given by the tree- level diagram with the Cf interaction shown in Fig. |l|(a). The result is 
simply 

zT 1(a) = -iCfk.k'. (6) 
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The second term in the low-momentum expansion is suppressed by k 2 /M^ compared to 
the leading order. It is given by the tree-level diagram with the Cf interaction shown in 
Fig. 0(b): 

iT m = -iC p A k 2 k • k' . (7) 

At order 12irk 5 /mM^, we have the one-loop diagram with two Cf interactions shown in 
Fig. |](c). The contribution of this diagram is 



iT m = {-iClf J 



d A q i q • k i q • k' 



where the dq$ integral was performed via contour integration. The remaining integral must 
be proportional to 5{j since no other vectors are available. Adding and subtracting k 2 in the 
numerator, we find 



k 2 — ie 



zT 1(c) = (Cf) 2 |^k-k' |y dqq 2 + k 2 J dq + k 4 J dq-^ 

= (Cf)^k-k' [L 3 + k 2 L 1 + ^ey (9) 



where L 3 and L\ are infinite constants. These two ultraviolet divergent terms can be ab- 
sorbed by redefining the low-energy constants C\ and Cf, respectively, which are already 
present in Ti( a \ and Tin,)- No new parameter enters at this order. The series proceeds in an 
obvious way. 

We can now match to Eq. (H) to relate the renormalized coefficients to the effective-range 
parameters. We find from Eq. @ that C\ = 127ra 1 /m, and from Eq. ([?]) that C\ = C\r\a\j2. 
After renormalization, T^) reproduces the third term in the low-momentum expansion of 
Eq. (|). Note that dia gram 0(c) cannot be renormalized by Cf alone even though it does 
not contain a C\ vertex. This observation has important consequences in the unnatural case 
with fine-tuning. 



B. Unnatural Case 



Now we turn to the more interesting case with a shallow p-wave state. In Refs. [§, |10 



it was shown that for a shallow s-wave state the leading-order contact interaction Co has 
to be treated nonperturbatively. In this case, Co is enhanced by a factor Mhi/Mi Q over the 
expectation Co ~ An/mMhi from naive dimensional analysis. Adding a new rung in the 
ladder forming the amplitude means adding an intermediate state (~ mk/Air) and a Co 
(~ An/mMio). Since the physics of the bound state is determined by k ~ M; Q , C has to be 
summed to all orders. 

For p waves matters are slightly more complicated. We have seen above that the renor- 
malization of the one-loop diagram with two C\ interactions requires tree-level counterterms 
corresponding to both the leading C\ and subleading C\ interaction. As consequence, at 
least the C\ and C\ interactions have to be treated nonperturbatively if a shallow p-wave 
state is present. 
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Bound (virtual) states are associated with poles in the S'-matrix on the upper (lower) 
half of the complex momentum plane. The characteristic momentum 7 of the bound/virtual 
state is given by position of the pole, \k\ = 7. For a shallow p-wave state with 7 ~ Mi Q , the 
magnitude of both the effective range and the scattering length must be set by M/ Q . This 
is a consequence of the renormalization argument from the previous subsection. Either C\ 
and C\ are both enhanced or they are both natural. Assuming that the higher terms in the 
effective range expansion are natural, the order of magnitude of the first three terms in the 
expansion is 

k 3 cot 5, ~ Mf Q + M lo k 2 + -r^—k 4 + . . . , (10) 

and the effective-range parameters scale as 

i~AC ^~M,„, and £ ~ jL . (11) 

Both C\ and C\ are enhanced over the expectation from naive dimensional analysis and 
scale as 

Cl~^T and Ci~J^-. (12) 

2 mMf 4 mMf D y ' 

Consequently, for momenta of order M\ Q neither interaction can be treated perturbatively. 
The shape parameter Vi, however, is of order I/Mm and its contribution is suppressed by 
Mi /M hi compared to the leading order. 

In the following, we will demonstrate that treating the Cf and C£ interactions to all 
orders is indeed sufficient for proper renormalization and, moreover, required to reproduce 
the physics of the shallow p-wave state. We will also work out the leading-order description 
of a shallow p-wave state. 

For convenience, we will not use the Lagrangian (^) but follow Ref. |2IJ and introduce an 
auxiliary field (the dimeron) for the two-particle state. The corresponding Lagrangian is, 



L = i/j 



id + 



, )/; JV> + Vi4[i9 + — -A 1 jd l + ^ (4(VViV) + H.c.) + . . . , (13) 

where the sign rjx = ±1 and the parameters g\ and Ai will be fixed from matching. This 
Lagrangian contains exactly the same number of parameters as the original Lagrangian (|j). 
Up to higher order terms, Eq. (|13|) is equivalent to Eq. (|j), as can be seen by performing 
the Gaussian path integral over dj. 

The bare dimeron propagator is given by 

tf>?(Pb, Pk = J? 6 " a ■ (14) 

Po ~ P /4m — Ai + it 

Summing the C\ and Cf interactions to all orders in the theory without the dimeron cor- 
responds to dressing the bare dimeron propagator with particle bubbles to all orders. This 
summation is shown diagrammatically in Fig. [| The full dimeron propagator is most easily 
calculated by first computing the self-energy from the particle bubble — i£i which up to 
overall factors is given by the one- loop diagram from Fig. |l|(c). We have 

d 4 q qiqj 



z£i(p ,p)ij = g I J 



(2TlY (PO _L a (P/ 2 +q) 2 i (PO n (p/2-q) 2 , iA 

t6lJ lt {l L3 + l L1 ( mpo ~ p2/4 ) + 1 ( mpo ~ p2/4 ) 3/2 } ' (15) 
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FIG. 2: The full dimeron propagator (thick shaded line) is obtained by dressing the bare dimeron 
propagator (double solid line) with particle bubbles (solid lines) to all orders. 



where L3 and L\ are infinite constants as in Eq. (^|). The full dimeron propagator now 
simply follows from the geometric series 

zDi(po,p) = i£>?(po,p)+*-D?(Po,p)HS 1 (p ,p))i£)?(p ,p)+... 

= iD° 1 { Po ,p)(l-i: 1 (p ,p)D 1 (p ,p)y 1 , (16) 

where the vector indices have been suppressed. Using Eqs. (|H[ fL5|) , we find 



■ / \ . s 12vr /12ttAi 12tt , 2 . n 

«A(Po,P)ij = -»dy — 2 2 l m ^o - P /4 (17) 

myf \r]imgi r)\m 2 gx v 7 

2 2 / 9 \ / 2 \ 'V^\ ^ 

-L 3 Li [mpo - p /4) - i [mp - p A) 

71 71 v 7 v y 



. 12tt / 12vrAf 12tt / 2 . x / 2/ x3/2\ 1 

where the last line defines the renormalized parameters Af and . 

The p-wave scattering amplitude is obtained by attaching external particles lines to the 
full dimeron propagator. In the center-of-mass system, (po,p) — (k 2 /m, 0), this leads to 



m y m(gi) z m z {gi* z 



[27r . k 2 cos9 f_l + r _i k 2 _ ik A 1 , (is) 



m \ ai 



from which the matching conditions can be read off easily. We see that, as advertised, two 
coefficients are necessary and sufficient to remove any significant cutoff dependence. 



C. Pole Structure 

In this subsection, we discuss the pole structure of the S-matrix in the unnatural case. 
Neglecting terms suppressed by Mi a /Mhi, the equation determining the poles is, from the 
amplitude (|IBD, 

- — + —k 2 - in 3 = 0. (19) 
di 2 

For defmiteness, we concentrate on the case ai,r± < that is relevant to na scattering. 
Other cases can be examined as easily. The solutions are one pole K\ on the positive 
imaginary axis and two complex-conjugated poles in the lower half-plane. They have the 
structure 

«i = ill and k± = £(7 ± £7) , (20) 
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FIG. 3: The pole structure of the S'-matrix for a p-wave resonance. 



where 



7i = | (l^l + 

a/3 / la 
7 = ~U H 



ai I 1 / 3 |rj 2 


+ hl 1/3 ) 


V 




ai 1/3 ri 2 


V 


2w 


2|a 1 | 1 /3 


|l/3| n |2 


" ) 







I \ V3 

108+ |ai||r 1 | 3 + 108^/1 + |a 1 ||r 1 | 3 /54j . (21) 

This pole structure is illustrated in Fig. |3|. This general structure remains qualitatively 
unchanged in the limit \r\\ — > 0. 

The p-wave contribution to the S-matrix can be written as 

c = e 2iSi = k + K i k + K + k + K- = _fc+j7i E-Ep- %T(E) 
1 6 " k — Ki k — k + k — k_ k - i'7i E - E + |r(.E) ' 

where we have defined 

k 2 t 2 + 'y 2 /"fiT 

E = -, E = ^L and r(£) = -W- , (23) 
2/x 2// V 2/i 

with /i the reduced mass of the system. The phase shift can therefore be written as 
Here 



5 s (S) = iarctanf^|j (25) 
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is the contribution from the bound state with binding energy B = jf/2fi. It changes by tt/2 
as the energy varies across B. 5 S (E) is a relatively smooth function of the energy E. The 
two complex-conjugated poles k± generate the resonance that is given by the second term 
in Eq. (|24]). This term changes by ir as the energy varies across E . 

In the case of s waves, the EFT determines in leading order the position of a shallow real 
or virtual bound state. In the p waves the physics is richer: the two leading-order parameters 
provide the position and width of a resonance (in addition to the position of a bound state). 

III. APPLICATION TO ELASTIC na SCATTERING 

We are now in position to extend the EFT for shallow p-wave states from the previous 
section to the n- 4 He system, including the spin of the nucleon. We calculate the leading- 
and next-to- leading-order contributions to low-energy elastic na scattering. First, we briefly 
review the structure of the cross section and scattering amplitude. 

A. Cross Section and Scattering Amplitude 

The differential cross section for elastic na scattering in the center-of-mass frame can be 
written as 

^ = \F(k,6)\ 2 + \G(k,6)\\ (26) 

where k and 6 are the magnitude of the momentum and the scattering angle, respectively. 
The so-called spin-no-flip and spin-flip amplitudes F and G can be expanded in partial waves 
as 

F(k, 9) = ^[(l + Vft+W+lf^P^cose), (27) 

l>0 

G(k,9) = J2[fi+(k)-fi-(k)} Pi (cos 9), (28) 
i>i 

where Pi is a Legendre polynomial and 

Pl{x) = {l-x 2 fl 2 j-P l {x). (29) 
The partial wave amplitudes fi± are related to the phase shifts 6i± via 



1 



2i8± 



1 



fl±{k) = m e 1 = ^ — (30) 



k cot Si± — ik 



The total cross section can be obtained from the optical theorem, 

a T = -r ImF(fc,0). (31) 
The T-matrix calculated in EFT is related to the amplitudes F and G via 

T = — [F + itr ■ hG) , (32) 
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where // = m a rriN / (m a + m^) is the reduced mass, n = k x k'/|k x k'| with k and k' the 
initial and final momenta in the center-of-mass frame, and cr = (<7i, 02, 03) is a three- vector 
of the usual Pauli matrices. 

For na scattering at low energies only the s and p waves are important. There is one 
s wave: Z± = 0+ with lj = S1/2, and two p waves: Z± = 1+ and 1— corresponding to 
h = P3/2 and P1/2, respectively. In the remainder of the paper, we use the Z± notation for 



the partial waves. In Ref. a phase-shift analysis including the 0+, 1—, and 1+ partial 
waves was performed and the effective-range parameters were extracted. The effective range 
expansion for a partial wave with orbital angular momentum I was given in Eq. (fj). The 
effective-range parameters extracted in Ref. fl9|] are listed in Table Q. The 1+ partial wave 
has a large scattering length and somewhat small effective range, as expected from Eq. (|i~l~l). 
Indeed, the phase shift in this wave has a resonance corresponding to a shallow p-wave state 
As a consequence, the 1+ partial wave has to be treated nonperturbatively using the 
formalism for shallow p-wave states developed in the previous section. In the 0+ wave, on 
the other hand, the scattering length and effective range are clearly of natural size. The 
0+ partial wave can be treated in perturbation theory. The situation is less clear in the 
1— wave. Although the pattern is similar to the 1+ wave, the phase shifts in the 0+ and 
1— partial waves show no resonant behavior at low energies fl9|j . We therefore expect that 
perturbation theory can be applied to the 1— partial wave as well. 

These points can be made slightly more precise. We can estimate the scales M\ and 
from the effective-range parameters. Using the parameters for the 1+ partial wave from 
Table [I], we find for M/ Q 50 MeV from the scattering length and 90 MeV from the effective 
range. The average value is M; G « 70 MeV. From the shape parameter, we extract ~ 260 
MeV. This is consistent with the hierarchy M lo -C M hi ~ ~ ^m^E a , where E a = 20.21 
MeV is the excitation energy of the a core PH, and suggests that our power counting is 



appropriate for the 1+ partial wave. We would expect that the scale of all effective-range 
parameters in the remaining channels is set by M^. Extracting the numbers, however, we 
find for M hi the scales 80 MeV from a 0+ , 280 MeV from r 0+ , 80 MeV from ai_, and 40 MeV 
from ri_. While some spread is not surprising given the qualitative nature of the argument, 
these numbers suggest that, even though the 1— phase shift is small, this partial wave might 
also be dominated by M[ . For the moment we will assume this is not the case and treat 
the 1— wave in perturbation theory. We can certainly improve convergence by resumming 
1— contributions. We return to this point in Sect. |J fD]. 



Partial wave l± 


a l± [fm 1+2i ] 


ri± [fm 1 2l ] 


V l± [fm 3 - 2 '] 


0+ 


2.4641(37) 


1.385(41) 




1- 


-13.821(68) 


-0.419(16) 




1+ 


-62.951(3) 


-0.8819(11) 


-3.002(62) 



TABLE I: The values of the scattering length ai±, the effective range r/±, and the shape parameter 
Vi± in elastic na scattering for the 0+, 1—, and 1+ partial waves from Ref. [19|. The numbers 
in parenthesis indicate the error in the last quoted digits. All values are given in units of the 
appropriate powers of fm as determined by the orbital angular momentum I of the partial wave. 
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B. Scattering Amplitude in the EFT 



A real test of the power counting comes only by calculating the amplitude at various 
orders and comparing the results among themselves and with data. In the following, we 
will compute na scattering to next-to-leading order in the EFT. For characteristic momenta 
k ~ Mi , the leading-order contribution to the T-matrix is of order I2ir /mMi . The EFT 
expansion is in M^jM^ and the NLO and N 2 LO contributions are suppressed by powers 
of Mi /Mhi and M^/M^, respectively. The parameters in the effective Lagrangian will be 
determined from matching to effective-range parameters. We then compare our results with 
the phase-shift analysis |19j and also directly with low-energy data. 

We represent the nucleon and the 4 He core by a spinor/isospinor N field and a 
scalar /isoscalar field, respectively. We also introduce isospinor dimeron fields that can 
be thought of as bare fields for the various Na channels. In the following we will employ s, 
d, and t, which are spinor, spinor and four-spinor fields associated with the Si/ 2 , P1/2, and 
p 3 /2 channels, respectively. 

The parity- and time-reversal-invariant Lagrangians for LO and NLO are 4 



C 



LO 



NLO 



id + 



2m r 



+ A f 



id + 



+ ^±Lt S t. N^(f)-(^N)(f) 



2m N 
+ H.c 



Ac. 



9o+ 



2(m a + m N 
t f S f • ^(N<P) + H.c. 

^ 2 



id + 



A 



1+ 



2(m a + m N ) 



t 



(33) 
(34) 



where r = (m a — m N )/{m a + m N ). The notation is analogous to that in Eq. (|13|). The 
S^s are the 2x4 spin-transition matrices connecting states with total angular momentum 
j — 1/2 and j = 3/2. They satisfy the relations 



SiS 



Sin 



^ijk^k 



5*j Sj 



4- 1 e f /2 



(35) 



where the jf^ 2 are the generators of the J = 3/2 representation of the rotation group, with 



T 3/2 j3/2 



7-3/2 



(36) 



These Lagrangians generate contributions in the 1+ and 0+ partial waves, 
contributions in N 2 LO, and the 1— partial wave enters first at N 3 LO. 
The propagator for the field is 



There are no 



iS^(p ,p) 



Po - P 2 /2m Q + it 



(37) 



4 We make a particular choice of fields here. The 5-matrix is independent of this choice. One can, for 
example, redefine the t field so as to remove the g' 1+ term. In this case, its contribution (see Eq. (|47 
below) is reproduced by a t* N<f> (+ H.c.) interaction with three derivatives. 
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while the nucleoli propagator is 



z£W(po,p w ' 



a/3 



Po - P 2 /2m N + it ' 



(38) 



In Eq. (|38"D, a and (3 (a and 6) are the incoming and outgoing spin (isospin) indices of the 
nucleon, respectively. The bare propagator for the 1+ dimeron is 



Po - p 2 /2(m a + m N ) - A i+ + ie 



(39) 



with a and (3 (a and b) the incoming and outgoing spin (isospin) indices of the dimeron, 
respectively. Note that 5 a p is a 4 x 4 unit matrix, since the dimeron carries j = 3/2. The 
bare propagator for the 0+ is slightly different because its kinetic terms do not appear until 
higher order: 



"WW- t- 

^04 



(40) 



with 5 a p now a 2 x 2 unit matrix. The bare propagator for the 1— dimeron is the same as 
for the 0+ dimeron, with the index 0+ replaced by 1— where appropriate. 

The leading contribution to the not scattering amplitude for k ~ M\ is of order Yl-n jmM\ 
and comes solely from the 1+ partial wave with the scattering- length and effective-range 
terms included to all orders. The next-to- leading order correction is suppressed by M^/M^ 
and fully perturbative. It consists of the correction from the shape parameter V\ + in the 1+ 
partial wave and the tree-level contribution of the scattering length ao+ in the 0+ partial 
wave. The 1— partial wave still vanishes at next-to- leading order. 

First, we calculate the leading-order T- matrix element T L °. As demonstrated for spinless 
fermions in the previous section, this is most easily achieved by first calculating the full 
dimeron propagator for the 1+ dimeron and attaching the external particle lines in the end. 
Apart from the spin/isospin algebra, the calculation is equivalent to the one for spinless 
fermions that was discussed in detail in the previous section. The proper self energy is given 
by 



iE 1+ (po,p)S 



9l + 



d d l (1 - rp/2) i (5j) /3T (l - r V /2) 3 {S 3 ) ia 5 a 



(27r) d ( 



f + l 



-igl + (S i Sj) / 3 a Sba 



_ (p/2+l)2 
1 1 2m a 

d d -H 



+ ie 



f PO / 



(P/2-1) 2 

2mjv 



+ ie 



(1 - rp/2Ml - rp/2), 



(2vr) a! - 1 po - (p 2 /4 + l 2 -rp- l)/2/i + ie 



(41) 



where we have performed the dlo integral via contour integration. Evaluating the remaining 
integral using dimensional regularization with minimal subtraction for simplicity, we obtain 



Si+(p ,p; 



a/3 



-5 a0 5, 



P<~>ab- 



6n 



2/i -Po + 



P 



-i 3/2 



2(m a + m N ) 



ie 



(42) 



Using Eq. (p!6|), the full dimeron propagator is then given by 



iD 1+ (p ,p 



)% = iVi + S aP S ab ipo - 2{m * +mN) ~ Ai + 



6tt 



-Po + 



P' 



2{m a + m N ) 



ie 



3/2 



ze . (43) 
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The leading-order T-matrix element in the center-of-mass system is obtained by setting 
(Po,p) = (k 2 /2fi, 0) and attaching external particles lines to the full dimeron propagator. 
This leads to 



T LO (k,cos9) = —k 2 (2cos6 + t(T -nsinfl) ( m+ ^±-Tj 1+ -^--ik 3 ) ' . (44) 



Using Eqs. (0) and ( |2"T| ) to (|52|), we find the matching conditions 

6ttA 1+ /i 2 #f + 

which determine the parameters A 1+ , and the sign rji + in terms of the effective-range 
parameters a\ + and r 1+ . Then, 

F^(k,9) = 2PcOS ^ 



-l/oi+ + r 1+ k 2 /2 - ik 3 ' 

— l/a 1+ + r 1+ k A /2 — ik 6 

At next-to-leading order, we include all contributions that are suppressed by Mi Q /Mhi 
compared to the leading order. These contributions come from the shape parameter V\+ 
and the s-wave scattering length cto+. Using the Lagrangian (^) and (p4|), we find for the 
T-matrix element: 

t nlo = Vo+9j+ ^g[+ 2?r k 6 (2 cos 9 + icr ■ n sin 9) 

A 0+ + 2^g\ + ii (-l/a l+ + r l+ k 2 /2-ik 3 f 1 ] 



The first term in Eq. ( f4T|) corresponds to the scattering length in the 0+ wave, while the 
second term corresponds to the 1+ amplitude with the shape parameter treated as a per- 
turbation. Using Eqs. (|) and (0) to ([32]), we find 

a 0+ = ~ o a and ^1+ = -3-2 ' 48 

Note that to this order i] 0+ , g 0+ , and A + are not independent and only the combination 
appearing in Eq. (|48D is determined. The next-to-leading-order pieces of F and G are then 



iNlo / 1 m _ , 2£; b cos0 



F WLO (M) = -a 0+ + 



4 (-l/ai+ + r 1+ P/2-ifc 3 ) 



2 



tNLO/K m _ Pi+ fc 6 sinfl 



C. Phase Shifts and Cross Sections in the EFT 

In order to see how good our expansion is, we need to fix our parameters. In principle we 
could determine the parameters by matching our EFT to the underlying EFT whose degrees 
of freedom are nucleons (and possibly pions and delta isobars) , but no core. Unfortunately, 
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FIG. 4: The phase shift for na scattering in the 1+ partial wave as a function of the neutron kinetic 
energy in the a rest frame. The dashed (solid) line shows the EFT result at LO (NLO). The filled 
circles show the phase-shift analysis |19| which the EFT at NLO reproduces exactly. The dotted 
line shows the contribution of the scattering length alone. 

calculations with the latter EFT have not yet reached systems of five nucleons 0. For the 
time being, we need to determine the parameters from data. For simplicity, we use the 
effective-range parameters from Table p] together with Eqs. fl4"5U4"8 ). 



In Fig. ||, we show the phase shifts for elastic na scattering in the 1+ partial wave as a 
function of the neutron kinetic energy in the a rest frame. The filled circles show the phase- 



shift analysis of Ref. | I5fl . The dashed line shows the EFT result at leading order. The 
LO result already shows a good agreement with the full phase-shift analysis. As expected, 
the agreement deteriorates with energy. NLO corrections improve the agreement: the EFT 
result at NLO shown by the solid line reproduces the phase-shift analysis exactly. If better 
data were available and a more complete phase-shift analysis were performed, some small 
discrepancies would survive, to be remedied by higher orders. 

The sharp rise in the 1+ phase shift past it/ 2 denotes the presence of a resonance. To 



LO, the pole structure of the S'-matrix is given in Sect. [II C| . We find 71 = 99 MeV, 7 = — 6 
MeV, and 7 = 34 MeV. Using Eq. (^3|), the position and width of the resonance are E = 0.8 
MeV and T(E ) = 0.6 MeV, respectively. The two virtual states that produce the resonance 
are indeed at \k\ ~ M\ . The real bound state, for reasons that cannot be understood from 
the EFT itself, turns out numerically to be at considerably higher momentum, where the 
EFT can no longer be trusted. This is consistent with the known absence of a real 5 He 
bound state. 

We also illustrate in Fig. [| an important aspect of the power counting. The dotted 
line shows the result from iterating Cf alone. In other words, it is the contribution of 
the scattering length only. This curve, which would come from a naive application of the 
power counting for s waves H |K]], does not correspond to any order in the power counting 
developed here, and clearly fails to describe the resonance near E n = 1 MeV. 
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FIG. 5: The phase shift for na scattering in the 0+ partial wave as a function of the neutron 
kinetic energy in the a rest frame. The solid line shows the EFT result at NLO; at LO this phase 
shift is zero. The filled circles show the result of the phase-shift analysis . 



In Fig. |j| we show the phase shifts for elastic na scattering in the 0+ partial wave as a 
function of the neutron kinetic energy in the a rest frame. In LO the phase shift is zero. 
The solid line shows the EFT result at next-to-leading order. The NLO result already shows 
good agreement with the full phase-shift analysis |19fl , depicted by the filled circles. 

The phase shifts in the 1— and all other partial waves are identically zero to NLO. The 
first non-zero contribution appears at N 3 LO in the 1— channel. All other waves appear at 
even higher orders. That they are indeed very small one can conclude from their absence in 
the phase-shift analysis fl9|j . 

Obviously, not all partial waves are treated equally in our power counting. In order to 
further assess if the power counting is appropriate, we compare the EFT predictions directly 
to some observables. In Fig. we compare the EFT predictions with data for the total 
cross section as a function of the neutron kinetic energy in the a rest frame. The diamonds 
are "evaluated data points" from Ref. In order to have an idea of the error bars from 
individual experiments we also show data from Ref. [M as the black squares. The dashed 



line shows the EFT result at LO which already gives a fair description of the resonance 
region but underestimates the cross section at threshold. The NLO result given by the solid 
curve gives a good description of the cross section from threshold up to energies of about 4 
MeV. 

We can also calculate other observables. As another example, we show in Fig. |7| the center- 
of-mass differential cross section at a momentum kcM — 49.6 MeV. (This corresponds to a 
neutron kinetic energy E n = 2.05 MeV in the a rest frame.) The diamonds are evaluated 
data from Ref. p2l 5 . The dashed line show the EFT results at LO, which is pure p wave. At 



In order to obtain the differential cross section from the NNDC neutron emission spectra we divide by 2tt 
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FIG. 6: The total cross section for na scattering in barns as a function of the neutron kinetic 
energy in the a rest frame. The diamonds are evaluated data from Ref. , and the black squares 
are experimental data from Ref. |23|]. The dashed and solid lines show the EFT results at LO 
and NLO, respectively. The dash-dotted line shows the LO result in the modified power counting 
where the 1— partial wave is promoted to leading order. 



NLO, shown as a solid line, interference with the s-wave term gives essentially the correct 
shape. 

If we carry out the EFT to a sufficiently high order, we will have included all terms used 
in the phase-shift analysis fli~9| , and more. At this order, the high quality of our fit is purely 
a consequence of the high quality of that fit. Note, however, that this is by no means true 
at the lower orders explicitly displayed above. In particular, it is perhaps surprising that 
our 1— wave does not appear until relatively high order. The fact that the EFT converges 
fast to data shows that the power counting developed here is reasonable. The 1— wave is 
further discussed in the next section. 



D. Further Discussion of the Power Counting 

As we have shown, the EFT describes the data pretty well at least up to E n = 4 MeV 
or so. One way to improve the convergence at higher energies is to take the scale of non- 
perturbative phenomena in the 1— wave as a low scale. We can modify the power counting 
and count the 1— parameters the same as the 1+ parameters. The LO Lagrangian from 



and multiply by the total cross section. 
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FIG. 7: The differential cross section for na scattering in the center-of-mass frame in barns/sr 
as a function of the scattering angle 9cm at a momentum kcM = 49.6 MeV. The diamonds are 
evaluated data from Ref. [22|. The dashed and solid lines show the EFT results at LO and NLO, 



respectively. The dash-dotted line shows the LO result in the modified power counting where the 
1— partial wave is promoted to leading order. 



Eq. ( |33D then has an additional term 



ALO 



9i 



id + 



2(m a + m N ) 
N^<f> - (^N)(f) 



d 

H.c. 



dV • \f(iV0) + H.c. 



(50) 



The calculation of the T-matrix for the 1— partial wave proceeds exactly as for the 1 + 
partial wave. The amplitudes F and G acquire the following additional contributions at 
leading order 



F ALO (k,9) 
G ALO (M) 



k 2 cos 6 



-l/ai_ +r^k 2 /2-ik 3 ' 

k 2 sin 6 
-l/ai_ +n_fc 2 /2 -ik 3 ' 



(51) 



In Fig. H, we show the phase shifts for na scattering in the 1— partial wave obtained in 
this alternative power counting. The filled circles show the result of the partial-wave analysis 



of Ref. |L9[ which is exactly reproduced by the leading-order EFT with the modified power 
counting given by the dash-dotted line. The dotted line shows the contribution of the 
scattering length only. The next-to-leading order in the modified counting cannot be easily 
computed at present because V\- is not known. 

The cross sections corresponding to the leading order in the modified power counting 
are shown by the dash-dotted curve in Figs. || and [7[ In the total cross section, promoting 
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FIG. 8: The phase shift for na scattering in the 1— partial wave as a function of the neutron 
kinetic energy in the a rest frame. The filled circles show the result of the phase-shift analysis [19] 
which is exactly reproduced by leading-order EFT with the modified power counting given by the 
dash-dotted line. The dotted line shows the contribution of the scattering length only. 



the 1— partial wave to leading order gives almost no improvement compared to the original 
counting except at higher energies, but even there the NLO result in the original counting 
gives better results. In the differential cross section at kcM = 49.6 MeV (corresponding to 
E n = 2.05 MeV in the a rest frame), the alternative power counting gives no improvement 
over the LO result (compare the dashed and dash-dotted lines in Fig. |7|). We did not find a 
significant improvement in the differential cross section over the leading order by promoting 
the 1— partial wave for neutron energies up to E n « 4 MeV. For reproducing the differential 
cross section, the interference between Sx/2 and p%/2 waves is much more important than the 
additional pi/% contribution. As a consequence, we deem the original power counting most 
appropriate for elastic na scattering at k ~ M\ . 

Finally, note that for k <C M\ the power counting has to discriminate between momentum 
k and the low-energy scale M; D . The p waves, for example, die faster than the s waves. That 
is the reason our results for the cross section in this region are not good until we get to 
NLO. It is easy to adapt the power counting for k -C M; G : in fact, the full amplitude — all 
waves, that is — can be treated in perturbation theory, as in Sect. |I A[ For more details, 
see Ref. M. 



IV. CONCLUSION AND OUTLOOK 

In this paper we have examined the problem of the interaction between a neutron and an 
a particle at low energies. We showed that a power counting can be formulated that leads 
to consistent renormalization. In leading order, two interactions have to be fully iterated. 
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These two interactions generate a shallow £>3/2-wave resonance near the observed energy and 
width. In subleading orders the phase shifts in all waves can be systematically improved. 
Observables calculated directly are very well reproduced. 

The crucial ingredient for the applicability of the EFT to bound states and resonances 
of halo type is their low characteristic energies. In this sense, the deuteron can be thought 
of as the simplest halo nucleus whose core is a nucleon. not scattering plays an analogous 
role here as np plays in the nucleons-only EFT. It is clear now how to extend the EFT 
to more complicated cores: one simply introduces an appropriate field for the core under 
consideration, extends the power counting to the relevant channels, and determines the 
strength of interactions order-by-order from data. 

With the parameters of the nucleon-core interaction fixed in lowest orders, we can proceed 
to more-body halos. The simplest example is 6 He. In addition to the na interaction, the nn 
interaction has also been determined from data. 6 He, like the triton, can be described as a 
three-body system of a core and two neutrons. The role of a three-body interaction can be 
addressed by renormalization group techniques [13|, 



Note that the EFT approach is by no means restricted to neutron halos. The Coulomb 



interaction can be included in the same way as in the nucleons-only sector [pi], allowing for 



the analysis of nuclei such as 8 B. Radiative capture on halo nuclei, such as p + 7 Be — *> 8 B + 7, 
can then be calculated much like n + p — > d + 7. 

Our approach is not unrelated to traditional single-particle models. In the latter, the 
nucleon-core interaction is frequently parametrized by a simple potential with central and 
spin-orbit components The parameters of the potential are adjusted to reproduce 

whatever information is accessible experimentally. In the EFT, we make the equivalent to 
a multipole expansion of the underlying interaction. The spin-orbit splitting, in particular, 
results from the different parameters of the dimeron fields with different spins. In the EFT 
the nucleon-nucleon interaction is treated in the same way as the nucleon-core interaction, 
mutatis mutandis. Contact NN interactions have in fact already been used in the study of 
Borromean halos |2Efl . It was found that density dependence, representing three-body effects, 
needed to be added in order to reproduce results from more sophisticated parametrizations 
of the NN interaction. In the EFT, the need for an explicit three-body force can be decided 
on the basis of the renormalization group before experiment is confronted. A zero-range 



model with purely s-wave NN and nucleon-core interactions was examined in Ref. p7 . 

The EFT unifies single-particle approaches in a model-independent framework, with the 
added power counting that allows for an a priori estimate of errors. It also casts halo nuclei 
within the same framework now used to describe few-nucleon systems consistently with 
QCD 0, ||. Therefore, the EFT with a core can in principle be matched to the underlying, 
nucleons-only EFT. Nuclei near the drip lines open an exciting new field for the application 
of EFT ideas. It remains to be seen, however, whether these developments will prove to be 
a significant improvement over more traditional approaches. 
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